arXiv:1505.07813v2 [cond-mat.mtrl-sci] 22 Dec 2015 


Large directional anisotropy in multiferroic ferroborate 

A. M. Kuzmenko/ V. Dziom/ A. Shuvaev/ Anna Pimenov/ M. Schiebl/ A. A. 

Mukhin/ V. Yu. Ivanov/ I. A. Gudim/ L. N. Bezmaternykh/ and A. Pimenov^ 

^Prokhorov General Physics Institute, 

Russian Academy of Sciences, 119991 Moscow, Russia 
^Institute of Solid State Physics, Vienna University of Technology, lOfO Vienna, Austria 
^L. V. Kirensky Institute of Physics Siberian Branch of RAS, 660036 Krasnoyarsk, Russia 

(Dated: December 23, 2015) 

Abstract 

One of the most fascinating and counter-intuitive recent effects in multiferroics is the directional 
anisotropy, the asymmetry of light propagation with respect to the direction of propagation. In 
such case the absorption in a material can be different for opposite directions. Beside absorption, 
different velocities of light for different directions of propagation may be also expected, which is 
termed directional birefringence. In this work, we demonstrate large directional anisotropy in a 
multiferroic samarium ferroborate. The effect is observed for linear polarization of light in the range 
of millimeter-wavelengths, and it survives down to low frequencies. The dispersion and absorption 
close to the electromagnon resonance can be controlled by external magnetic held and is fully 
suppressed in one direction. By changing the geometry of the external held, samarium ferroborate 
shows giant optical activity, which makes this material to a universal tool for optical control: 
with a magnetic held as an external parameter it allows switching between two functionalities: 
polarization rotation and directional anisotropy. 
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I. INTRODUCTION 


Multiferroics are materials which exhibit electric and magnetic order simnltaneouslyi^-. 
Due to the coupling of electric and magnetic effects, these materials show a strong potential 
to control electricity and magnetism and, more generally, the properties and propagation 
of light. An unusual way to influence the light beam in a material is provided by direc¬ 
tional anisotropy. This effect is the inequivalence of forward and backward directions of 
propagation and it may be divided into directional dichroism and directional birefringence. 
Directional dichroism results in an asymmetry of the absorption coefficient. Directional 
birefringence may be dehned as different velocity of light for forward/backward direction, 
or, equivalently, by difference in the refractive index. 

One way to obtain the asymmetric transmission is a material is given by the magnetochiral 
anisotropy^. This effect is physically close to the Faraday rotation^ and results in non¬ 
equivalent transmission for light propagating parallel or antiparallel to the magnetisation of 
the sample. Experimental demonstrations of magnetochiral anisotropy exist in a series of 
systems^^- including ferromagnets^ and multiferroicsl^^— . 

Further possibility to break the symmetry with respect to propagation direction may be 
realized in multiferroics and with magnetization perpendicular to the light propagation. In 
multiferroics with electric polarization P perpendicular to the magnetization M a nonzero 
toroidal momentii^ appears, T = P x M, which may be parallel or antiparallel to the 
propagation direction. Directional anisotropy based on toroidal moment has been recently 
demonstrated in multiferroics!^'— and is realised in present work. 

More strictly, directional anisotropy may be dehned in terms of non-reciprocity and is 
equivalent to an asymmetry with respect to interchanging the source and detector. Mathe¬ 
matically, non-reciprocal effects in optics can by quantihed via a property called reactioi>!^ 
which measures the difference between backward and forward geometries. In several specihc 
cases, the reciprocity (i.e. equivalence of forward and backward directions) can be proven 
rigorously, e.g. for scalar helds!^ or for media with symmetric susceptibilities and without 
magnetoelectric terms22,. A recent review about the topic may be found in Ref. 2]|. Ear¬ 
lier results on non-reciprocal optics^Si'— have been obtained in a classical magnetoelectric 
material Cr203. 
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The propagation of light in the matter is governed by the material equations: 


B — /io/iH + 

D = yTo/fox®’”!! + £o^E. 

Here B is the magnetic flux density, D is the electric displacement held; E and H are 
the electric and magnetic helds; jl and Xme/em are the matrices of electric permittivity, 
magnetic permeability and magnetoelectric susceptibilities, respectively. As may be shown 
by rigorous calculation^'^ the necessary conditions for reciprocal propagation are given by 

e = e'^, fl = fF, Xme = -Xem ' (2) 

Here ()^ denotes the transposed matrix. 

One well-known case of a non-reciprocal transmission is given by the Faraday rotation. 
Here the effect is realised by antisymmetric off-diagonal elements of permittivity or perme¬ 
ability. As will be shown later, in case of samarium ferroborate the non-reciprocity is due 
to symmetric off-diagonal terms in the magnetoelectric susceptibility (xTy)- 

An instructive question within the scope of the present work is: What is the symmetry 
operation exchanging the source and detector, i.e. which operation reverses the direction of 
light? If we consider plane electromagnetic waves as a standard experimental tool, both time 
inversion (t —t) and space inversion (r —r) would reverse the propagation direction. 
In case when linearly polarized waves are the eigenwaves of the problem (as in experiments 
below), both time and space inversions exchange the direction but the linear polarization 
of the wave is preserved. The situation gets signihcantly different if circular (or elliptically) 
polarized waves are the propagating eigenmodes in the sample. We recall that besides 
changing the propagation direction, space and time inversions differently influence the helds 
of the electromagnetic wave: time inversion preserves the direction of the electric held, but 
inverts the magnetic held. The space inversion does just the opposite: the electric held is 
inverted and the magnetic held is preserved. This diherence has fundamental consequence for 
the circular polarization. As may be easily shown, the space inversion interchange clockwise 
and counterclockwise polarizations (i.e. the handedness), but the time inversion does not. 
From the dehnition of reciprocity it can be shown that the rotation sense of circular waves 
must be kept after inverting the propagation direction. Therefore, in practical experiments 
with inverted propagation of light, space inversion symmetry does not seem to be adequate 
at least in case of circular (or elliptical) polarizations. 
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The propagation of light within a medium is characterized by a wave vector and - for 
a given propagation direction - may be reduced to a complex refractive index. The real 
part of the refractive index is responsible for the collection of a phase shift during the 
propagation. In other words, it reflects the ratio between the light velocity in the media and 
in vacuum. The imaginary part of the refractive index is responsible for the energy loss and it 
is termed absorption coefficient. Directional anisotropy of the absorption coefficient is called 
directional dichroism. Strong directional dichroism in multiferroics arises due to coupling 
of electric and magnetic order— In addition, specific symmetry conditional^ should 
be fulhlled in such experiments, e.g. the orthogonal mutual orientation of magnetization, 
electric polarization, and propagation direction. 

Electromagnetic waves are solutions of Maxwell equations, where six components of elec¬ 
tric and magnetic held and the wave vector (propagation vector) are unknowns. Using two 
Maxwell equations, VD = 0 and VB = 0, two held components can be removed from the 
system. Here D and B are electric displacement and magnetic hux density, respectively. 
Using boundary conditions, the problem reduces to a four-dimensional linear problem for 
eigenvalues and eigenvectors. Four possible eigenvalues represent four possible propaga¬ 
tion constants in the media. In general case, all eigenvalues are diherent. However, due 
to inversion symmetry in the majority of the problems, two pairs of solutions are equal in 
absolute values and opposite in sign. They correspond to two possible polarizations of the 
electromagnetic wave in the sample with two (forward and backward) propagation directions 
for each wave. Non-reciprocal ehects arise when the eigenvectors are diherent for opposite 
directions. 

Here we present a material that can be considered as a universal tool for the control of the 
optical properties. Using solely an external magnetic held as a parameter two diherent ehects 
may be switched and modihed: (i) giant optical activity^! (polarization rotation) and (ii) 
giant directional birefringence and dichroism for linearly polarized light. For intermediate 
orientations both ehects are mixed realizing a general case of four-colored transmissionr^. 

Magnetoelectric rare-earth ferroborates^*^ RFe 3 (B 03)4 (R = rare earth ion) represent a 
newly discovered material class with strong magnetoelectric coupling. Especially in ferrobo- 
rates with R = Sm, Ho colossal magnetic held-induced changes in the dielectric constant have 
been observed^i^. In the case of Smo, 5 Lao. 5 Fe 3 (B 03)4 it has been proven experimentally 
that an intrinsic magnetoelectric excitation are responsible for the observed ehects^'^'^. 
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FIG. 1: Crystal structure of Smo^^Loo^^Fe^fBOsJi and directional birefringence, a - Basic struc¬ 
tural elements in Smo. 5 Lao. 5 Fe 3 (B 03 ) 4 . The direction of Fe-spins is given in the geometry with 
external magnetic field parallel to the 6-axis, b - Definition of directional birefringence. The 
backward beam (red) is only weakly refracted by the sample. 


Such excitations in magnetoelectric materials are called electromagnons^*^, and they are 
dehned as magnetic excitations that interact with the electric component of electromagnetic 
radiation. 

II. EXPERIMENTAL DETAILS 

Spectroscopic experiments in the terahertz frequency range (40 GHz < v < 1000 GHz) 
have been carried out in a Mach-Zehnder interferometer arrangement^ which allows mea¬ 
surements of the amplitude and the phase shift in a geometry with controlled polarization 
of radiation. Theoretical transmittance curves^S for various geometries were calculated from 
the susceptibilities within the Berreman formalisirt^. The experiments in external magnetic 
helds up to 7 T have been performed in a superconducting split-coil magnet with polypropy¬ 
lene windows. Static dielectric measurements have been done using commercial impedance 
analyzer equipped with a superconducting magnet. 

Large single crystals of Smo, 5 Lao, 5 Fe 3 (B 03 ) 4 , with typical dimensions of ~ 1 cm, have 
been grown by crystallization from the melt on seed crystals. 
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III. RESULTS 


A. Samarium Ferroborate 


The material used in the present study is Smo.5Lao,5Fe3(B03)4. This compound is closely 
similar to the pure samarium ferroborate^^i^, SmFe3(B03)4. Doping with lanthanum has an 
advantage of suppressing the growth of domains with different symmetry. The presence of 
domains in the sample would reduce the value of the magnetoelectric non-reciprocal effects. 
The crystallographic structure of Smo.5Lao,5Fe3(B03)4 is shown in Fig. la. It contains 
two interacting subsystems given by Sm^^/La^’*' and Fe^^ ions. The iron subsystem orders 
antiferromagnetically below Tn = 34 K with an easy-plane magnetic structure oriented 
perpendicularly to the trigonal c-axis. Although the Sm^"*" moments play an important role 
in the magnetoelectric properties of Smo,5Lao.5Fe3(B03)4, they probably do not order up to 
the lowest temperatures. 

Static electric polarization in multiferroic ferroborates can be explained by symmetry 
arguments and by taking into account that Fe^"*" moments are oriented antiferromagnetically 
within the crystallographic a 6 -plane^>^. Within the topic of the present work the terms 
governing the ferroelectric polarization along the a and 6 -axis (or x and j/-axes) are of basic 
importance. For the R 32 space group of borates these terms are given by 




Li-Li 


Py ~ ‘2LxLy . 


( 3 ) 


Here L = Mi — M2 is the antiferromagnetic vector with Mi and M2 being the magnetic 
moments of two (antiferro-)magnetic Fe^"*" sublattices. Further details of the symmetry anal¬ 
ysis of the static mag netoelectric effects in Smo.5Lao.5Fe3(B03)4 can be found in Appendix 
and in Refs. 39 - 41 |. 

The simple expression Eq. ([ 3 ]) allows to understand the behavior of static and dynamic 
properties in external magnetic fields. For magnetic fields along the a(a:)-axis an antiparallel 
orientation of electric polarization with respect to the a-axis is stabilized. In this case 
Smo.5Lao,5Fe3(B03)4 reveal strong electromagnon mode^Ii^ which can be excited separately 
via electric, e\\b{y), or magnetic, h\\a{x), channels. Here h and e are the ac-components of 
the electromagnetic radiation. As has been shown recently—, in the geometry B||a samarium 
ferroborate exhibits strong optical activity, i.e. the polarization of the incident radiation can 
be rotated on a controllable way. 


6 





B. Non-reciprocal transmission 


The geometry with magnetic field along the b{y)-axis is very promising from the point 
of view of non-reciprocal effects. External field B ||6 stabilizes the magnetic configuration 
with the magnetic moments oriented along the a-axis {Ly = 0 , ^ 0 ). In agreement 

with Eq. ([ 3 ]), the static polarization is oriented parallel to the a-axis. The electromagnon is 
excited simultaneously via the electric and magnetic channel: h\\a{x), e\\b{y). In this case the 
magnetoelectric coupling starts to distinguish between two possible propagation directions. 
This effect arises because electric polarization, P||a, is oriented perpendicular to the induced 
magnetization M||B|| 6 . Consequently, Smo,5Lao,5Fe3(B03)4 gets a toroidal^ moment, T = 
P X M, which may be parallel or antiparallel to the propagation direction^. The toroidal 
moment allows the existence of strong directional birefringence in Smo,5Lao.5Fe3(B03)4. 

The susceptibility matrices in Eq. ([ 1 ]) can be written in a 

simplified form (see Appendix): 
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As shown in more details in the Appendix, for an electromagnetic wave propagating 
along the z-direction the complex refractive indexes of the forward {w^) and backward (n“) 
solutions are different: 


n 
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\/ Sy^x i 


a 


xy 


( 5 ) 


Here axy = x 7 y “ xTyX 7 z/f^z, Sy = Sy + {x^yf j^z, and fix = ^x + (x^JV/^^ are the 
renormalized material parameters; Sy = 1 + Xyy is the dielectric permittivity along the 
?/-axis, Hx = ^ + xTx and = ^ + xS are the magnetic permeabilities along the x-axis 
and x-axis, respectively. In present approximation The solutions in Eq. ([ 5 ]) are 

written for a linearly polarized electromagnetic wave with a polarization e||?/, h||x. The 
solution for the perpendicular polarization with e||x, h\\y is trivial and can be written as 
= y/Ex^y Importantly, two possible solutions for the propagating wave are linearly 
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FIG. 2 : Non-reciprocal transmission in SmQ^^Lao^^FesfBOs)^. Left panels; non-reciprocity is 
achieved by the inversion of magnetic field, a - Forward/backward transmission, b - Non-reciprocal 
phase shift, <plv. Right panels: non-reciprocity is achieved by explicit exchange of source and detec¬ 
tor. c - Forward transmission for two opposite directions of external magnetic field, d - Backward 
transmission. The geometry of the experiment is given in the inset. Symbols - experiment, solid 
lines - fits according to the Berreman model as described in the Appendix. 

polarized. Therefore, no polarization rotation is expected for the geometry with B ||6 and 
the two linear polarizations do not mix. This is in contrast to a related geometry with B||a 
revealing optical activity^. 

Typical transmittance spectra in the geometry with non-reciprocal effects (B|| 6 ) in the 
frequency range of our spectrometer are shown in Fig. O Switching between forward and 
backward geometry is experimentally achieved by inverting the direction of the external 
magnetic held, i.e. the direction of the induced magnetization, M||B. We recall that in 
agreement with Eq. (| 3 ]) the direction of electric polarization does not change, but the toroidal 
moment M x P does. Therefore, the inversion of the external magnetic held is equivalent 
to the inversion of the propagation direction. From the point of view of symmetry, two 
experimentally relevant cases are: M x P 'I'j' k and M x pnk. 








In order to check explicitly, that the exchange of the sonrce and detector does lead to 
the same results, the corresponding experiments have been carried out. The results of these 
experiments are presented in the right panels of Fig. [21 Only the absolute values of the 
transmission have been measured in this case as the Mach-Zehnder interferometer— can not 
be easily inverted. Fig. | 2 t shows the transmittance in the geometry with forward propagation 
direction and with two opposite values of external magnetic held. We state that similarly to 
the spectra in the left panels, the electromagnon is strong for one direction of the magnetic 
held only. The transmittance spectra for the backward direction are presented in Fig. [2)d. 
Here the situation is just the opposite to that in the panel c: the held values with excited 
and silent electromagnon interchange. 

We note that no Faraday rotation is expected in present case, because the Voigt geometry 
with B||fe, k||c has been used. 

Figs. |2^,c,d demonstrate the dependence of the transmission intensity on the propaga¬ 
tion direction, which is equivalent to a directional dichroism. In Smo.5Lao.5Fe3(B03)4 the 
directional dichroism is observed only close to the electromagnon resonance frequency, 
z/q ~ 146 GHz (H = 6 T). Already for frequencies ±10 GHz apart from the resonance, 
forward and backward transmission intensities coincide within experimental accuracy. 

C. Directional birefringence 

Figure 2b shows the diherence in phase shift (optical thickness) in 
Smo.5Lao.5Fe3(B03)4 for forward/backward propagation. Again, close to the electro¬ 
magnon frequency strong diherence in optical thickness is observed between two directions. 
An important contrast to the asymmetric transmission in Fig. 2a is a broadband character 
of the observed ehect, i.e. substantial variation in optical thickness is detected in the full 
frequency range of the present experiment. 

Figure [ 3 | demonstrates that two solutions for the propagating waves are linear and do not 
mix. In the geometry with h||a clear electromagnon excitation is observed close to ±6 Tesla, 
as in this conhguration electric and magnetic excitation channels are active. We note here 
again the contrast between positive and negative directions of the magnetic held. In a simple 
interpretation using Eq. (jS]) for the propagation constants, one can say that for positive helds 
the magnetoelectric contribution is added to the ’’conventional” electromagnon mode. For 
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FIG. 3 : Absence of polarization rotation. Magnetic field dependence of the transmittance spectra in 
Smo.sLao.sFes (603)4. Positive/negative values of magnetic field is equivalent to forward/backward 
direction of propagation, respectively, a - Amplitude of the transmitted signal, b - Relative 
changes of the optical thickness. Open symbols (blue and black) - magnetoelectric geometry of 
the experiment with both, electric and magnetic channels, excited. Full symbols - silent geometry: 
the electromagnon is not excited. Solid lines - fits as described in the Appendix. Absence of the 
electromagnon excitation in the silent geometry (B||/i|| 5 ) demonstrate zero polarization rotation. 

negative magnetic fields the magnetoelectric susceptibility is subtracted from the refractive 
index. The characteristic values of the susceptibilities are close to typical numbers, for which 
the electromagnon is nearly suppressed in the backward direction. 

The silent geometry in Fig. [ 3 ] with B||h||6, e||a shows no excitation in the spectra. This 
observation supports the simple model presented above, in which no electromagnon and no 
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FIG. 4 : Directional birefringence and dichroism in SmQ,^LaQ,^Fe-i(BO'i)ii. a - Refractive index in 
Smo.5Lao.5Fe3(603)4 for forward (blue) and backward (red) propagation of the millimeter-wave 
radiation, b - Directional birefringence, c - Absorption coefficient in Smo.5Lao.5Fe3(603)4 for 
forward and backward propagation, d - Directional dichroism. Symbols - experiment, solid line - 
fits according to the model described in the text. 

magnetoelectric contribution is present. Most importantly, this demonstrates the absence 
of the polarization rotation in the conhguration B||6, k||c. That is, active and silent polar¬ 
izations do not mix and remain linear. The observed effect remains correct as long as the 
spectra are dominated by a low-frequency electromagnon. The effect of the high-frequency 
electromagnoii^'^ at about 320 GHz can be neglected. As may be expected, if the interac¬ 
tion of both electromagnons is substantial, certain effects of the polarization rotation would 
appear. 

IV. DISCUSSION 

Figure m summarizes the central results of this work, i.e. directional birefringence and di¬ 
rectional dichroism in Smo.5Lao.5Fe3(B03)4. As described in more details in the Appendix, 
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the refractive index {n) and the absorption coefficient (k) have been obtained from the 
transmission and phase shift spectra shown in Fig. [2l As the formulas do include the inter¬ 
ferences within the sample surfaces, Fabry-Perot-like oscillations seen in the transmittance 
spectra (Fig. | 2 ]) are substantially suppressed in the data of Fig. SJ As expected from the 
measured spectra of the optical thickness, the refractive index in Smo.5Lao.5Fe3(B03)4 is 
strongly direction-dependent. Moreover, a substantial difference in n is present at all fre¬ 
quencies shown in Fig. 01 The difference is also clearly seen in Fig. 0j3 which directly presents 
the directional birefringence (i.e. An = — n~). Because the electromagnon frequency 

strongly depends upon an external magnetic held, vq ~ B, the birefringence can be varied 
by magnetic held and even changes sign on crossing the resonance frequency. 

In agreement with the data in Fig. 0 ]b, the directional dichroism remains nonzero even 
at static frequencies. According to Eq. (jS]) the static value is 

An(0) = n+(0) - n-(0) ^ 2x^;(0) . (6) 

We note that the static value of the refractive index is ill-dehned with increasing wavelength 
and Eq. ([6]) is probably relevant down to microwaves only. 

In a contrast to the directional birefringence, the directional dichroism in 
Smo,5Lao.5Fe3(B03)4 is limited in frequency (see Fig. 0 b, d). Although the dichroism is 
strong close to the electromagnon frequency, is gets unmeasurably small ±10 GHz apart 
from the resonance frequency. 

A remarkable result seen in Fig. 0 ]is the nearly full suppression of the dispersion, n(z/), 
and absorption, k{u), for the backward propagation direction. In view of Eq. ([ 5 ]) for this 
set of parameters the magnetoelectric susceptibility almost completely suppresses the elec¬ 
tromagnon contribution. To analyze this behavior in some more detail we modify Eq. ([ 5 ]) 
taking into account the frequency dependencies of susceptibilities as given in the Appendix. 
We further assume in the hrst approximation that Ae\LF{oj)\ -C Soo, Afi\LF{oj)\ <C 1 , 
Afi -C Ae, and carry out the calculations to the linear order of Lf{uj) only. Here Ae 
and A/i are electric and magnetic contributions of the electromagnon, respectively; Lf{oj) 
is the Lorentzian function and £oo is the high-frequency dielectric permittivity. The hnal 
approximate expression may be written as: 
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FIG. 5: Model calculations of directional birefringence, a - Refractive index for forward propagation 
direction for typical range of the parameters in Smo.sLao.sFea(603)4. a - Same for backward 
propagation. Note the factor 8 difference in scale between a and b. 

From Eq. (( 71 ) we see that the effect of the magnetoelectric susceptibility on the electro- 
magnon is equivalent to a weighted sum or difference of electric and magnetic contributions. 
Therefore, we expect complete suppression of the electromagnon contribution if the condi¬ 
tion 




( 8 ) 


is fulhlled. As Ae: is large for small magnetic helds (Ae ~ 30 ) and is suppressed quadratically 
in high helds Ae ~ l/B"^ j^, some optimum value of the electromagnon suppression is 
expected in moderate helds, as observed in Figs. Ej 01 In such cases the directional anisotropy 
in Smo.5Lao,5Fe3(B03)4 is the strongest. 

In order to prove the estimates in Eq. ([8]) by rigorous calculations, we plot in Fig. [ 5 ] 
the model curves for the refractive index in Smo,5Lao.5Fe3(B03)4 for forward and backward 
propagation directions, respectively. Typical absolute values of the characteristic parameters 
used in Fig. [ 5 ] have been taken from the hts to the data above and from the results of 
previous experiments^Ii^. From the scale diherence in Fig. [ 5 ] we see that the amplitude 
of the Lorentzian is suppressed by at least a factor of An^/An~ ~ 8 in the full range of 
parameters. In agreement with estimates by Eq. ([8]) the suppression of the electromagnon 
for backward direction is the strongest for magnetic helds in the range i? ~ 5 — 8 T, as 
observed in the experiment. 

Compared especially with hrst experiments^^— on non-reciprocal transmission in Cr203 
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the effect in samarium ferroborate is extremely strong and may approach absolute asymme¬ 
try in absorption. Such values are due to strong coupling of electric and magnetic order and 
is a characteristic feature of multiferroic materials with strong electromagnonii^— . We note 
that in the data in Fig. Il^,b no electromagnon mode is observed for backward polarization 
within experimental accuracy. 


V. CONCLUSIONS 

Directional birefringence and dichroism have been investigated in doped samarium fer¬ 
roborate at millimeter-wave frequencies. We demonstrate strong directional anisotropy close 
to the resonance frequency of electromagnon excitation. The strength and position of the 
relevant electromagnon can be varied in external magnetic field. The non-reciprocity is 
demonstrated using two alternative approaches: i) by the inversion of the external magnetic 
field and ii) by explicit exchange of the source and detector. Nearly full suppression of the 
electromagnon in one propagation direction is observed. Approximate range of parameters 
to obtain strong directional anisotropy is investigated. 
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Appendix 

Symmetry analysis of the susceptibilities 

Pure and La-substituted SmFe 3 (B 03)4 belong to R32 space group with easy-plain anti¬ 
ferromagnetic structure at low temperatures. At low temperatures and in external mag¬ 
netic fields along the y-axis the following electric and magnetic configuration is stabi- 
lized2^>2^‘22,^— M||j/, L||a;, P||x. Here M ~ /xqH is the induced magnetization, L is 
the antiferromagnetic vector, and P is the static ferroelectric polarization. In this case the 
magnetic structure corresponds to 2^ point group. From the symmetry arguments^ the 
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following elements of the magnetoelectric susceptibility matrix are allowed: X^y’yT,zx,xz- W^e 
note that in the dynamic case several additional terms (see below) may become nonzero. 

Therefore, in order to obtain the relevant elements of the susceptibility matrices, a model 
of the statio^>^>^‘^ and dynamio 2 ^> 2 S>^ properties of Smo.sLao.sFea( 603)4 has been utilized. 
The magnetoelectric energy $me(M, L,P) relevant for the present analysis can be written 
ag^di; 


d)^e(M,L,P) 


—Cl [PxLyLz — PyLxLz] — C 2 \Px {L^. — Ty) — 2PyLxLyj + ... . (9) 


The dynamic susceptibilities in Smo.sLao.sFes( 603)4 are governed by two electromagnon 
modes^IiS^iiS which we denote as F-mode and AF-mode. In the geometry B\\b two modes are 
coupled and the susceptibility matrices ;\;'"(a;), x'"®(a'), x®'"(ci;), x®(a;) are obtained so 


Landau-Lifshitz equations for dynamic magnetic variables AM, AL similar to Ref. 27l |: 
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where the components Xxy,yx,xz,zx^ X7lyy,zy, ^%,yx are of purely dynamic origin and 

they are proportional to iuj. The absence of the components x7z,yz,zz, Xl7,zy,zz, Xlz,yz,zx,zy 
results from neglecting of a higher-order magnetoelectric term PzLzLx{L‘l. — SL^). 

The AF mode has been detected close to 320 GHz at low temperatures^. The resonance 
frequency of the F-mode is close to zero (~ 5GHz) but increases roughly linearly in external 
magnetic helds^^. In case which is relevant for the present experiment, the interaction 
between two electromagnon modes can be neglected. In this approximation and in the 
vicinity to the resonance frequency up the susceptibility matrices may be simplihed to the 
form given in Eq. (jl]) of the main text. 
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Data processing and transmission spectra 


We start from Maxwell equations which for plane waves (E, H ~ e**^^^* can be written 


as: 


X E = = ^n;[/io(i + , (11) 

-ik X H = iooD = ia;[y'£ohoX^™'H + £o(i + X^)E] . (12) 


Here 1 is the identity matrix. 

In general, Eqs. (11111121) can be reduced to eigenvalue and eigenvector problems giving 
the wavevectors and the amplitudes of the propagating waves. 

The electromagnetic held inside the sample has the form of a plane wave exp{i{kxX + 
kyU + kzZ — ojt)). The values of k^ and ky are conserved at the boundaries and they are 
determined by the geometry of the problem. The value of k^ depends on the properties of 
the sample and is obtained solving the Maxwell equations within the sample follo wing the 
Berreman method^. Further details of the calculations may be found in Refs. 27lJ36| . In 
present experiments the light propagates along the 2 ;-direction, i.e. the k vector may be 
written as k = (0, 0, k). 

The solutions of Maxwell equations, Eqs. (Illlll2j) . are given by two linearly polarized 
modes. One solution with h|||/, e||x is trivial and reveals no dependence on the propagation 
direction: 

k^\/2'K = ±y/e^ (13) 


Here = 1 + Xxx dielectric permittivity along the x-axis and = 1 + y™ is the 

magnetic permeability along the y-axis. 

The propagation constant for a perpendicular polarization with h||x, e||j/ does show non¬ 
reciprocal effects and it reads: 


kfX/2n = ttxy ± \/eyHx (14) 

Here a^y = x7y “ xTyX7z/f^z^ + ixTyf/f^z, and jlx = fix + (x™JVh^ are the 

renormalized material parameters; Sy = 1 + Xyy is the dielectric permittivity along the y- 
axis, Hx = ^ + xTx and = ^ + X^x are the magnetic permeabilities along the x-axis and 
; 2 -axis, respectively. In present approximation ~ hz- For a single-domain sample, the 
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relevant electrodynamic parameters in Eq. flT^ are obtained as^i; 


Sy — 1 + Xyy — ^OO + • Lpiuj) 

/Xj; = 1 + X™ = 1 + A/r ■ Lpioj) 

X7z = {iuj/uJo)^f^- Lf{uj) 

Xxy = P^/AfiAe- Lpiu) 

X7y = -i—^/ApAe ■ Lf{uj) 
ojf 

Here £00 is the high-freqnency contribntion to the permittivity, uj = 2 ttv is the angular 
frequency, uq, Ap, and Ae are the resonance frequency, magnetic, and dielectric contribution 
of the low-frequency electromagnon, respectively. The factor p{H) = H/^JIP^PiW^He 
reflects the changes of the magnetic structure with increasing magnetic field and it becomes 
unity in the helds exceeding 5-10 kOe. Here H'j^ and Hp are anisotropy and exchange helds, 
respectively. 

We utilize the Lorentzian form for the electromagnon resonance with 

Lf{u:) = [up - - 1 - iug) , 


with g being the resonance width. 

The normalized impedance of the medium for the solution given by Eq. (ITTl) . pq/eq = 
Ey/Hx, does not depend upon the propagation direction and may me written as: 


z = 



(15) 


The transmission amplitude in case of a plane parallel sample can be calculated explicitly 


as 


(1 - r^)e 


2\Ak^d 


(16) 


_ y,2^i{k++k-)d 

Here r = (z — l)/(z -|- 1) is the reflectivity on the sample surface. 

We utilized Eq. ffTOj) to obtain the £t curves in Figs. [2]in We note that in Eq. ffT6|l only 
the term is asymmetric with respect to the propagation direction. This means that 
relative birefringence and dichroism may be directly obtained from the ratio /t~ without 
knowing the details of the experiment. In order to rigorously invert the transmission data 
in Fig. [21 the values of r are missing in Eq. (IT6|) . However, as we can get from the fitting 
procedure, r{u) varies only by ±3% in the full range and by ±0.2% far from the resonance. 
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Therefore, and without loosing the accuracy, r{u}) can be chosen as a constant. In order to 
improve the approximation, r{u) has been taken from the hts to the transmission data in 
Fig. [2J With known r{uj) the inversion of Eq. fllbp and the calculation of k~^ and k~ from 
t~^ and t~ is straightforward. The results of this procedure are shown as symbols in Fig. HI 
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